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III-V growth and surface conditions strongly influence the physical structure and resulting optical 
properties of self-assembled quantum dots (QDs). Beyond the design of a desired active optical 
wavelength, the polarization response of QDs is of particular interest for optical communications 
and quantum information science. Previous theoretical studies based on a pure InAs QD model 
failed to reproduce experimentally observed polarization properties. In this work, multi-million 
atom simulations are performed to understand the correlation between chemical composition and 
polarization properties of QDs. A systematic analysis of QD structural parameters leads us to 
propose a two layer composition model, mimicking In segregation and In-Ga intermixing effects. 
This model, consistent with mostly accepted compositional findings, allows to accurately fit the 
experimental PL spectra. The detailed study of QD morphology parameters presented here serves 
as a tool for using growth dynamics to engineer the strain field inside and around the QD structures, 
allowing tuning of the polarization response. 
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Introduction: Semiconductor nanostructures are be- 
ing more and more applied to several optoelectronics 
technologies ranging from lasers^ to optical amplifiera^l 
or single photon sourceP^ where they have successfully 
overcome critical challenges such as extremely low thresh- 
old, high speed response, or entangled photon emission, 
respectively. In most of these applications a crucial 
parameter is the polarization response, typically mea- 
sured in terms o f th e degree of polarization (DOP = 
TE-TM/TE-l-TMJElllor the TM/TE raticP. Understand- 
ing how this feature is related to the quantum dot (QD) 
structural symmetry and composition could be very help- 
ful in its tuning. For example, semiconductor optical 
amplifiers (SOAs) for telecommunications require engi- 
neered QDs for isotropic polarization behavioi'^EEl. Po- 
larization control is also crucial for other applications, 
such as polarization-entangled photons emitted by single 
QDsS and polarization sensitive applications of vertical 
cavity surface emitting lasers^ ^. 

The structural and electronic properties of III-V 
QDs epitaxially formed by the self-ordering Stranski- 
Krastanov (SK) process are strongly affected by surface 
and growth conditions. Several techniques such as sur- 
face structuring^, nanostructure engineering via strain 
coupling^ -^'-'^^ ^^, and tuning of growth conditional^ are 
currently being studied for controlling QD distribution 
and geometric symmetry. As a common trend, epi- 
taxial InAs nanostructures exhibit asymmetry not only 
in the growth plane (with respect to [110] and [110] 
directions JpEI but also in the vertical plane (along [001])!^ 
due to their preferential flat, lens shape. The vertical con- 
finement leads to a strong compressive biaxial strain sup- 
pressing the light-hole (LH) component in the valence- 



band states leaving mostly the heavy hole (HH) compo- 
nent and resulting in predominantly in-plane (TE) po- 
larized emission from the interband transitiongSl. 

InAs QDs obtained via the SK process have been 
found to be significantly influenced by In-Ga intermix- 
ing and In-segregation effects du ring the capping and 
post-growth annealing processea^^tiS. Several composi- 
tion profiles have been proposed and different investiga- 
tion techniques employing high resolution transmission 
electron microscopy, x-ray diffraction, photoelectron mi- 
croscopy or scanning probe microscopy have been used 
to exactly map InAs nanostructuregi^HSil As recently re- 
viewed by Biasiol and HeurP^l^ results in the literature do 
not lead to an unified model, where the actual composi- 
tion profile of the QDs is strongly related to the growth 
conditions. Their review highlights a common tendency 
that the chemical composition of a typical SK QD has 
gradients both along the growth and the in-plane direc- 
tions: the In composition increasing from base to top due 
to In segregation effects and decreasing from the center 
towards the edges in the lateral directions due to In-Ga 
intermixing effects. However, such a complex structure 
has not been considered so far to theoretically understand 
the polarization respo nse of InAs QDs. 

Previous theoreticaPESMI studies of QD polarization 
response are based on a pure InAs type QD composi- 
tion profile, thus significantly limiting their accuracy and 
leading to discrepancy between theory and experiment. 
A clear example of such a discrepancy is the failure to 
reproduce the large values of the e xperimentally mea- 
sured DOP by using both k • jJlHESl and atomistic tight 
binding methods^l. Thus a better understanding of the 
correlation between the QD structure and the degree of 
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FIG. 1: (a) Atomic force microscopy (AFM) image of a single QD layer sample before GaAs capping revealing a dot density 
of 3.5x10^° dots/cm'^ with average height of 5nm as indicated in the inset above the figure, (b) Normalized PL intensity is 
plotted for TEioo (broken lines) and TMooi modes (solid lines), from experimental measurements and theoretical calculations 
as a function of wavelength. The calculations are done for a pure InAs and an Ino.7Gao.3As random alloy quantum dot. The 
comparison of the experimental measurement and the theoretical calculations clearly indicates a significant difference for the 
TM/TE ratios, (c) Plots of the ground state wavelengths (GSW) and TM/TE ratio as a function of the In-composition of the 
QD. As the In composition decreases from 100% (pure InAs QD) to 55% (an alloyed Ino.55Gao.45As QD), the GSW decreases 
drastically whereas the TM/TE ratio only slowly increases. The experimental values of the GSW and TM/TE ratio are marked 
as dashed lines. Clearly a single In-composition in the uniform composition model is unable to simultaneously reproduce 
experimental values of the GSW and TM/TE ratio. 



polarization remains an outstanding challenge. 

In this letter, the polarization response of InAs QDs 
is theoretically studied by atomistic simulations, intro- 
ducing a compositional model capable of fitting experi- 
mental measurements both of the electronic transitions 
and of the TM/TE ratio. The actual complex composi- 
tion and geometry of SK InAs QDs is mimicked by a two 
composition QD model, reproducing the experimentally 
measured polarization behavior of a single QD layer and 
highlighting the relevance of atomic scale processes like 
segregation and intermixing. 

Experimental procedure: The QD samples used 
in this study were grown on semi-insulating GaAs sub- 
strates, by a COMPACT 21- Riber Molecular Beam Epi- 
taxy (MBE) system equipped with a reflection high en- 
ergy electron diffraction (RHEED) gun to monitor the 
surface evolution in-situ during growth. After growth of 
a GaAs buffer layer at 600° C, the substrate temperature 
was lowered down to 500° C and QDs were formed by cov- 
ering the buffer with 2.8 MLs of InAs. The 2D-3D growth 
mode transition is demonstrated by the RHEED pattern 
evolving from streaky-like to spot-like after deposition 
of 1.7 MLs of InAs. Afterwards, dots were immediately 
capped by a GaAs spacer layer grown at the same low 
temperature. The single layer QD sample was then cov- 
ered with a 20 nm GaAs cap terminating the structure. 
An uncapped single QD layer sample was also grown un- 
der the same conditions and its morphology analysed by 
atomic force microscopy (AFM), providing a dot den- 
sity of 3.5x10^° dots/cm^ with average height of 5nm, as 
shown in Figure llja). The GaAs capped QDs were es- 



timated by transmission electron microscopy (TEM) to 
have a dome-like shape with base diameter w 15nm and 
height sa 5nm (see inset above Fig. [l|^a)). 

For investigating the polarization behavior, samples 
were excited from the top with a cw Ar-|- laser (A — 514 
nm). The room temperature photoluminescence (PL) 
signal from the cleaved edge of the samples was first col- 
lected by a long focal lens (200 mm) and it was then fil- 
tered by a linear polarizer and focused by a second lens 
into the monochromator. 

Theoretical model: The th eoretical modeling is per- 
formed using NEMO S-EEIISI. NEMO 3-D is an atom- 
istic simulator based on the valence force field (VFF) 
methocpSl for strain calculations and the twenty band 
sp^d^s* tight binding modeP^ for the electronic struc- 
ture. NEMO 3-D is a multi-scale simulator capable of 
performing multi-million atom simulations for realistic 
QD dimensions surrounded by large GaAs buffers to 
model the long range impact of the strain and piezoelec- 
tric potential. This tool has already been used to model 
QD structures pr oviding r esults in good agreement with 
experimental datcPSEBEol Piezoelectric potentials, both 
linear and quadratic, are calculated by solving Poison's 
equation according to a published recipe^^l and are in- 
cluded in the calculations of the electronic spectra. The 
inter-band optical transition strengths are calculated us- 
ing Fermi's golden rule, with the squared magnitude of 
the optical matrix elements found by summing over the 
spin degenerate states. The polarization dependent TE 
and TM spectra are calculated along the [100] and [001] 
directions respectively, as a cumulative sum of optical 
transitions between the lowest conduction band energy 



Schematic of intermixed QD 
60nm 



(a) GaAs Capping Layer 

InxzGaAs 



IhxiG 

\ 




O.Snm 
InAsW 



^■nk \ 5nm 

> 



15nm 
GaAs Substrate 





(b) (B,H) 


0.40 


(10,5)o/^9'5) 


•B 0.35 
cc 


(9,4b^ 


S 0.30 


(10,4y^11,5) 


0.25 


j^^ Experiment 

(9,3)^11,4) 


0.20 


(10,3/ 



0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 
Aspect Ratio (H/B) 




11401 ' ' • ' ' • 

0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 



Aspect Ratio (H/B) 



FIG. 2: (a)Scheinatic diagram of the two composition QD model, placed on top of a 0.5nm thick InAs wetting layer (WL). 
The QD has overall dimensions of base=15nm and height=5nm in accordance with the TEM''''. It consists of two regions. The 
central In-rich region with Ina:2Gai_2,2As composition (x2 > 80%) is of size base=B and height=H. The outer In-poor region is 
made up of In^iGai-^^iAs (xl < 40%) material, (b) The TM/TE ratio as a function of the aspect ratio (AR=H/B) for various 
values, (B, H), in units of nm. The compositions xl and x2 are 40% and 100%, respectively, (c) The GSW as a function of 
the aspect ratio (AR=II/B) for various values of H and B. The compositions xl and x2 are 40% and 100%, respectively. The 
numbers in nm along the solid and dotted lines indicate the change in the value of GSW between the adjacent points along 
each line. 



level (El) and the highest four valence band energy levels 
(HI, H2, H3, and H4), where each transition strength is 
artificially broadened by multiplication with a Gaussian 
distribution centered at the wavelength of the transition—. 
The highest four hole energy levels are chosen here in- 
stead of only the top most level (HI), because the valence 
band states are closely spaced on the energy scale (Hl- 
H4 < kT « 26meV at T=300K) and multiple hole states 
contribute to the ground state optical emission peak at 
the room temperatureSl. 

Model 1: uniform compositions: We start our 

analysis by first assuming a uniform material profile 
throughout the QD geometry as in previous theo retical 
studies of the polarization response^ISlinilMiIll] Dq^^ 

are assumed to be lens-shaped with base diameter and 
height set at 15nm and 5nm, respectively, in accordance 
with the experimental TEM analysiP^I. Figurefnb) plots 
the experimental normalized PL spectra taken under the 
two polarization conditions (TE and TM) and compares 
them with the theoretically calculated spectra for both 
the ideal case of a pure InAs QD and an Ino.rGag.a As ran- 
dom alloy QD. The InGaAs alloy configuratio n has be en 
suggested by some previous theoretical studiea^^EIEII to 
mimic the In-Ga intermixing effect. 

The pure InAs QD model clearly fails to reproduce the 
experimentally measured spectra. The TM/TE ratio is 
calculated to be 0.097 which is significantly lower than 
the experimentally measured value of saO.265. The peak 
wavelength in the PL spectra, hereafter referred to as 
ground state wavelength (GSW), is also overestimated 
by Ril23nm. This is consistent with the earlier theo- 
retical studiepESMl where a pure-InAs-QD type model 
predicted TE-dominant optical emissions and failed to 



match the experimentally measured PL spectra. Lower- 
ing the average In composition inside the QD to 70% to 
mimic the effect of In-Ga intermixing during the capping 
process increases the TM/TE ratio to 0.115 due to the 
reduced biaxial strain, but at the same time blue shifts 
the GSW by 72nm with respect to the experimental value 
of 1226nm. To further investigate the impact of lowering 
the In composition of the QD, we gradually reduce its 
value from 100% (pure InAs QD) to 55% (Ino.55Gao.45As 
QD) and plot the corresponding values of the GSW and 
TM/TE ratio in Figure fl|c). As the In composition 
reduces, the GSW decreases to 1088nm, crossing the ex- 
perimental value (shown by the dotted line) for In ^84%. 
However, the TM/TE ratio increases very slowly, only 
reaching a value of 0.114 for In=55%. This clearly in- 
dicates that an uniform In-composition, cither as pure 
InAs or some InGaAs alloy composition, cannot be se- 
lected as the chemical composition of the QD when seek- 
ing agreement between theory and experiment for both 
the GSW and TM/TE ratio. Therefore, a more complex 
compositional model is required in order to understand 
the experimental results. 



Several ex perimenta P^ ' ^^ ' ^ -^J^ and 

theoreticaP^E^I^^IlIIlI] studies have shown that In- 
Ga intermixing and In-segregation effects significantly 
influence the chemical composition of QDs during the 
growth of the GaAs capping layer. Thus we hypothesize 
that a quantitative agreement between theory and 
experiment can be achieved if some inter-diffusion of Ga 
atoms inside the QD region is incl uded in the model. 
The most widely suggested profilcpJ^llMHlIl consists 
of an In-rich region along the vertical direction and 
a decreasing In-composition laterally from the center 
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FIG. 3: (a)Plots of the TM/TE ratio as 
a function of the compositions xl and x2 
as described in the QD schematic of Fig. 
[2ra). The plots are for two dimensions of 
the core region: (B, H) = (lOnm, 4nm) and 
(llnm, 4nm). The arrows indicate an in- 
crease/decrease of the TM/TE ratio with 
respect to a decrease in the compositions xl 
and x2. (b) Plots of the GSW as a function 
of the compositions xl and x2. The arrows 
indicate increase/decrease of the GSW with 
respect to a decrease in the compositions xl 
and x2. 



to the border. Such a complex configuration can be 
theoreticaUy approximated by defining two composition 
regions inside the QD: an In-rich core surrounded by a 
low In-concentration region. Following these considera- 
tions, we now analyse the polarization response of such 
a QD model and explore the structural parameters to 
reproduce the experimental measurements. 

Model 2: non-uniform composition: In order to 
get the most reliable configuration, as close as possi- 
ble to the experimental one^ a simplified double region 
scheme is adopted, consisting of an In-rich inner core of 
Inj,2Gai_2,2As (x2 typically > 80%) material, surrounded 
by an In-poor thin region of luajiGai-^^iAs (xl typically 
<40%), as schematically shown in Figure 2(a). These 
ranges for the compositions xl an d x2 h ave been chosen 
from previous experimental results^IESl. The actual size 
of the In-rich core, defined by inner region base width 
(B) and height (H) as shown in Figure l2ja), strongly 
depends on the growth dynamics of the capping layer. 
The overall size of QD obtained from TEM'^'^ places con- 
straints on the values of B and H: H < 5nm and B < 
15nm. On the basis of this QD model, we perform a sys- 
tematic investigation by changing the parameters B, H, 
xl, and x2 and comparing the corresponding GSW, TE, 
and TM modes with the experimental measurements. Fi- 
nally, it is worth noting that the two cases studied earlier 
in the uniform composition model, an InAs QD and an 
Ino.7Gao.3As QD, can be included by this two composi- 
tion model by assuming xl=x2=100% and xl=x2=70%, 
respectively. 

Variations of B and H: First, in order to investigate 
the impact of the parameters H and B on the polarization 
response, we choose x2=100% and xl=40%. Such an 
assumption is motivated by considering that the central 
core of InAs nano-islands is always reported as being In- 
rich and that the investigations reported so far indicate 
that the maximum value for the In-composition can reach 
1009p]. 

Figures ^h) and ^c) plot the TM/TE ratio and the 
GSW, respectively, as a function of the aspect ratio (AR 
— H/B) of the inner core. The target experimental values 
are shown as horizontal dotted lines. The TM/TE ratio 
in Figure ^b) linearly increases with AR when H is in- 
creased for a fixed value of B. This is because an increase 



in H will result in a reduction of the biaxial strain close 
to the center of the QD which will reduce the splitting 
between the HH and LH bands. As a result, the magni- 
tude of the TM-mode will increase, thus increasing the 
TM/TE ratio. However, our calculations show that the 
slope of the TM/TE ratio as a function of AR=II/B is 
different for different values of B. In general, it decreases 
as B increases, indicating a decreasing impact of H on the 
TM/TE ratio. This is consistent with a previous theo- 
retical studjBS where hole energy levels in smaller QDs 
were found to be more sensitive to changes in QD AR 
as compared to larger QDs. We want to highlight here 
that our theoretical calculations show that even changes 
of l-2nm in B or H values, which can actually be induced 
by highly controlled growth dynamics, result in a drastic 
change in the polarization properties of the QD sample. 

The plots of GSW as a function of AR in Figure W[c) 
exhibit a linear dependence of GSW on H for a fixed 
value of the B (shown by the solid lines), as well as a 
linear dependence of the GSW on the value of B for a 
fixed value of H (shown by the dotted lines) . The analysis 
of Figure [2jb) and (c) indicates that two sets of values 
for (B, H), namely (lOnm, 4nm) and (llnm, 4nm), give 
transition wavelengths and a TM/TE ratio close to the 
experimental values. A slight change in the sets of values 
of B and H ((lOnm, 5nm) and (9nm, 5nm)) leads to a 
GSW value close to the experimental one, but introduces 
a discrepancy in the corresponding TM/TE ratios. We 
therefore choose (B, H) to be (lOnm, 4nm) and (llnm, 
4nm) to analyse in the following sections the impact that 
varying xl and x2 has on the calculated values of the 
GSW and TM/TE ratio. 

Variation in xl : The outer shell composition xl was 
assumed above to be 40% in the study of the effect of the 
B and H parameters on the QD optical properties. How- 
ever, this composition can have a lower value depending 
upon the growth dynamics during the growth of the GaAs 
capping layer. We therefore decrease this composition to 
30% and 20% for both (lOnm, 4nm) and (llnm, 4nm) 
dimensions, while keeping x2 at 100%. Figures [sFa) 
and p[b) show the impact of such variations in xl on 
the TM/TE ratio and GSW, respectively. It is clear that 
decreasing the value of xl has opposite impact on the 
TM/TE ratio and GSW: the value of TM/TE ratio in- 



creases and the value of the GSW decreases. Following 
these trends for the dimensions (B, H) plotted in Figures 
WCb) and W[c), we deduce that a decreasing GSW with re- 
spect to xl will only bring (llnm, 5nm) and (lOnm, 5nm) 
closer to the experimental GSW. However, the TM/TE 
ratio for these two dimensions is already higher than the 
experimental value (see Figure ^b)), so a correspond- 
ing increase in the TM/TE ratio will further reduce their 
agreement with the experimental TM/TE ratio. There- 
fore, we deduce that decreasing xl below 40% reduces 
agreement with the experimental values (either GSW or 
TM/TE ratio) for all of the dimensions (B, H) and hence 
we require the poor In-composition region comprising QD 
outer shell to have a composition close to 40%. 

Variation in x2: Next, we analyse the impact of the 
last unknown parameter x2 in our proposed two composi- 
tional model (see Figure [sja)), keeping xl=40%. Figures 
Islb) and Islc) plots variations in the TM/TE ratio and 
GSW, respectively, when x2 is decreased from 100% to 
90% and 80%. Contrary to xl, a decrease in x2 has sim- 
ilar impact on both the TM/TE ratio and GSW: both 
decrease as the composition x2 decreases. From Figures 
Wih, c), these trends imply that the values of both the 
TM/TE ratio and GSW will become closer to the ex- 
perimental values only for the two dimensions, (llnm, 
5nm) and (lOnm, 5nm). However, the GSW for the di- 
mensions (lOnm, 5nm) is very close to the experimental 
GSW whereas the TM/TE ratio for this dimension is sig- 
nificantly different from the experimental value. Hence a 
collective agreement of both the GSW and TM/TE ratio 
with the experimental values is not possible for a decrease 
in x2. For example, at x2=90%, the GSW decreases to 
1187nm (w34nni below the experiment value) whereas 
the TM/TE ratio is stiU 0.31 (saO.045 above the experi- 
mental value). For the (B, H) = (llnm, 5nm) case, if x2 
decreases to 90%, the GSW and TM/TE ratio decrease 
to 1208nm and 0.27, respectively, which are in reasonable 
agreement with the corresponding experimental values of 
1221nm and 0.265. 

We complete our analysis of QD compositions by con- 
sidering simultaneous decrease in xl and x2 and finding 
that their cumulative impact on the GSW and TM/TE 
values is roughly equal to the sum of their individual ef- 
fects. We find that the TM/TE ratio is more sensitive 
to x2 as compared to xl, and therefore any simultaneous 
decrease in the compositions xl and x2 will reduce the 
values of both parameters (GSW and TM/TE ratio). 

QD structural parameters: From the preceding 
discussions on the effect of B, H, xl, and x2 in the QD 
structure of the Figure l2ja), we conclude that for our 
experimentally measured QD sample, the most effective 
model is represented by the following parameters: B w 
llnm, H « 4-5nm, xl « 40%, and x2 « 90-100 %. The 
calculated PL spectra for based on these sets of param- 
eters are compared with the experimental PL spectra in 
Figure |4] It is worth noting that the proposed struc- 
tures roughly reproduce the complex composition pro- 
file experimentally found by many authorP^l^ including 
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FIG. 4: Normalized PL spectra is shown from the experimen- 
tal measurement (black lines) and compared with the theo- 
retical calculations (green and red lines). Both TE (dotted 
lines) and TM (solid lines) components are plotted. A slight 
difference between the GSW of the calculated TM and TE 
modes is due to the fact that the TE mode is dominantly 
from El-Hl transition and the TM mode is dominantly from 
E1-II3 and E1-H4 transitions. Such a difference will be diffi- 
cult to observe in the experiment, because of the broadening 
of the PL spectra which are typically collected from multiple 
QDs. 



the anisotropic In-Ga intermixing behavior with respect 
to the crystallographic direction, with a larger in-plane 
inter-diffusion region (nearly 2nm around the base of the 
dot) and a negligible inter-diffusion at the top (0-lnm) 
where segregation predominated. 

Relaxed strain energies: The QD systems have 
been relaxed using the VF F mo del to reach a minimum 
strain energy configuratiorP^ESl. The comparison of the 
relaxed strain energies for the various QD systems under 
study may provide additional insight to find a more likely 
experimental geometry. Figure [Slplots the relaxed strain 
energies as a function of the total number of In-atoms for 
the various configurations after the VFF minimization is 
achieved. We consider two cases of the uniform com- 
position model: an Ina;Gai_a;As QD with decreasing In- 
composition (see Figure [lie)) and a pure InAs QD with 
decreasing base diameter. We also consider two cases 
of the non-uniform composition model: (B, H, xl=40%, 
x2=lG0%) and (B, H, xl=40%, x2=90%). Clearly, the 
relaxed strain energies belonging to the two layer conr- 
position configurations are lower than the strain energies 
of the uniform composition configurations for a similar 
number of In atoms. This suggests that a two composi- 
tion model is more favourable as a proposed experimen- 
tal geometry when compared to a uniform composition 
model, in accordance with our earlier findings. 

Strain profile analysis: Our model 2 agrees with 
the experimental TE and TM PL spectra. However a 
comparison of strain profiles is required to further under- 
stand the difference between the calculated PL spectra 
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FIG. 5: Plots of the total strain energy of the relaxed QD 
systems as a function of the total number of In atoms in the 
system. Both the uniform and non-uniform composition mod- 
els have been considered. The relaxed strain energies for the 
non-uniform composition configurations are smaller than the 
relaxed strain energies of the uniform composition configura- 
tions with similar number of In atoms. 



for models 1 and 2. In Figures p[a-c), we compare three 
QD configurations: (i) a pure InAs QD with base=I5nm 
and height=5nm, (ii) an Ino.7Gao.3As alloy QD of the 
same size, and (iii) a two-layer composition as in Fig- 
ure[2]^a) with B=llnm, H=4nm, xl=40%, and x2=100% 
which gives the best agreement with experiment. For ref- 
erence, the previously calculated values of the GSW and 
TM/TE ratio are also listed underneath. To understand 
the shifts in the values of the TM/TE ratio and GSW in 
(b) and (c) with respect to (a), we compare the hydro- 
static and biaxial strain plots along the [001] direction 
through the center of the QDs in Figures |6jd) and (e) , 
respectively. For the Ino.7Gao.3As alloy QD, the strain 
values inside the QD region are randomly distributed due 
to the disordered composition of the alloy. We mark the 
average values of the hydrostatic and biaxial strains by 
the dotted horizontal lines. 

The hydrostatic strain shifts the conduction and va- 
lence band edges, therefore increasing the band gap, 
whereas the biaxial strain is mainly responsible for con- 
trolling the splitting between the HH and LH valence 
band edges^. The strain driven band edge shifts in Ta- 
ble IT] are calculated from simple analytical expressions 
involving InAs/GaAs deformation potentials^Sl. 

Blue shift of the GSW: The blue shifts in the GSW 
for (b) and (c) with respect to (a) indicate increase of 
the optical gaps (Eg), which can be understood as a 
cumulative affect of the relaxed hydrostatic strains and 
the increase of the band gaps due to the reduced aver- 
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FIG. 6: (a-c) Schematics of three QD models under study: 
(a) A pure InAs QD, (b) A random alloy Ino.7Gao.3As QD, 
and (c) A QD with InAs core surrounded by an Ino.4Gao.6As 
region, (d, e) Plots of the Hydrostatic (e/f = e^x + ^yy + ^zz) 
and Biaxial [es = ^xx + ^yy — 2ezz) strain components along 
the [001] direction through the center of the QDs. The solid 
lines, open circles, and dotted lines are for the QD models 
shown in (a), (b) and (c), respectively. The horizontal dotted 
lines are marked to indicate the average values of the strain 
inside the QD region for the Ino.7Gao.3As alloy QD model. 



age In compositions of the QDs"^^. By comparing the 
hydrostatic strain profiles in Figure p^ (d), we deduce 
that the presence of InGaAs alloy, in general, relaxes 
the strain which can be attributed to reduced QD lattice 
mismatch with the GaAs buffer. As a result, the band 
gaps for (b) and (c) will be smaller when compared to 
(a) due to smaller strain driven shifts in the conduction 
and heavy hole band edges. For example, from Table [ij 
the strain driven shifts predict that the band gap for the 
Ino.7Gao.3As QD is 373.7-231.6=142. ImeV smaller than 
the InAs QD. However if the relative increase (w310meV) 
in the band gap of the Ino.7Gao.3As QD due to the re- 
duced average In composition is also added, the cumula- 
tive change is 310-142. 1=167. 9meV increase with respect 
to the InAs QD, consistent with the increase of the opti- 
cal gap. To summarize, the relaxation of the hydrostatic 
strain due to presence of the InGaAs alloy in (b) and (c) 
introduces reductions in the band gaps which are smaller 
than the corresponding band gap increase due to the re- 



TABLE I: The impact of the hydrostatic (en) and biaxial 
(en) strains on the conduction (Ec), heavy hole (Ehh), and 
light hole (Elh) band edges are given for the QD models 
shown in Figures|6[a, b, and c). The band edge deformations 
{SEcB, SEhh, and SF,lh) are calculated from the analyti- 
cal expressions given in the Ref. 29. S'Ega=S'EcB-HH and 
SF,h=SEhh-lh represent the changes in the band gap and 
the HH/LH splittings due to the strain. 
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(a) -0.0842 -0.1536 427.7 54 

(b) -0.06 -0.148 304.8 73.2 

(c) -0.0818 -0.1449 415.5 48.6 



-222.4 373.7 278.4 
-193.2 231.6 266.4 
-212.2 366.9 260.8 



duced average In compositions, and hence the cumulative 
impact of these two affects results in overall increase of 
the optical/band gaps (blue shifts of the GSWs) for (b) 
and (c) with respect to (a). 

Biaxial strain relaxation increases TM/TE Ra- 
tio: Figure J6[e) compares the biaxial strain components 
and Table ^provides the values of the corresponding 
changes in the HH and LH band edges for the three QD 
configurations. When compared to the pure InAs QD, 
the biaxial strain component is reduced for both the In- 
GaAs alloy QD and the 2-composition QD model. This 
reduces the splitting between the HH and LH bands, thus 
increasing the LH component in the top most valence 
band states. A past theoretical investigatiorr^ has also 
reported an increased HH/LH intermixing for an InGaAs 
QD when compared to a pure InAs QD. The increased 
LH character in turn enhances the TM componeniP, as 
we compute for configurations (b) and (c). The small- 



est HH/LH splitting responsible for the enhanced TM 
component occurs for the two-layer composition model. 

Conclusions: In conclusion, multi-million atom sim- 
ulations are performed to understand the correlation be- 
tween the chemical composition of self-assembled QDs 
and their electronic properties by reproducing the ex- 
perimentally measured polarization properties. A single 
composition QD model, being either a pure InAs QD 
or an alloyed InGaAs QD, failed to reproduce the ex- 
periment. To fully understand the experimental polar- 
ization behavior, the model must take In-Ga intermix- 
ing effects into account, by representing the QDs with 
variable compositions. Based on a systematic analysis 
of the dependence of the TM/TE ratio on the various 
QD morphology parameters, we propose a two composi- 
tion model that accurately reproduces the measured PL 
response. The model gives results consistent with the 
experimental compositional findings and highlights the 
strong anisotropy of atomic scale phenomena like inter- 
mixing and segregation affecting the polarization behav- 
ior of these nanostructures. These results could represent 
a tool for using growth dynamics to engineer the strain 
field inside and around the QD structures, allowing tun- 
ing of the polarization properties, a critical parameter for 
several challenging applications. 
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